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The statistical coalescence model for the production of open and hidden charm is considered within 
the canonical ensemble formulation. The data for the J/ip multiplicity in Pb+Pb collisions at 
158 A-GeV are used for the model prediction of the open charm yield. We find a strong enhancement 
of the open charm production, by a factor of about 2-4, over the standard hard-collision model 
extrapolation from nucleon-nucleon to nucleus-nucleus collisions. A possible mechanism of the open 
charm enhancement in A+A collisions at the SPS energies is proposed. 
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The charmonium states J /ip and tp' have been mea- 
sured in nucleus- nucleus (A+A) collisions at CERN SPS 
over the last 15 years by the NA38 and NA50 Collab- 
orations. This experimental program was motivated by 
a suggestion ]l| to use the J/tp as a probe of the state 
of matter created in the early stage of the collision. In 
this approach a significant suppression of J/tp produc- 
tion relative to Drell-Yan lepton pairs is predicted when 
going from peripheral to central Pb+Pb interactions at 
158 A-GeV. This is originally attributed to the formation 
of a quark-gluon plasma, but could be also explained in 
microscopic hadron models as secondary collision effects 
(see [El and references therein). 

The statistical approach, formulated in Ref. B, as- 
sumes that J/ip mesons are created at hadronization ac- 
cording to the available hadronic phase-space similar to 
other (lighter) hadrons. We call this a hadron gas (HG) 
model. The HG model offers a natural explanation of 
the proportionality of the J / tp and pion yields in nuclear 
collisions. 

Recently a picture of the J/ tp creation via coalescence 
(recombination) of charmed quarks was developed within 
several model approaches (4[j^] . Similar to the HG model 
H , the charmonium states are assumed to be formed at 
the hadronization stage. However, they are produced 
as a coalescence of created earlier c and c quarks and 
therefore the multiplicities of open and hidden charm 
hadrons are connected in the coalescence models. In 
Ref. H the charm quark-antiquark pairs are assumed 
to be created at the early stage of A+A collision and 
the average number of directly produced c-c pairs, N^ r , 
is fixed by the model consideration based on the hard 
scattering approach. This requires the introduction of a 
new parameter in the HG approach - the charm enhance- 
ment factor 7 C [Ej . This is analogous to the introduction 
of the strangeness suppression factor "f s || in the HG 
model, where the total strangeness observed is smaller 
than its thermal equilibrium value. Within this approach 
the open charm hadron yield is enhanced by a factor j c 
and the charmonium yield by a factor 7^ in comparison 
with the equilibrium HG predictions. 



As the total number of charmed hadrons is expected to 
be smaller than unity, even for the most central Pb+Pb 
collisions at the SPS energy, an exact charm conserva- 
tion within the canonical ensemble (c.e.) should be im- 
posed. This has been done in Ref. ||. Note also that 
the c.e. formulation was successfully used in Ref. Jl(J 
to calculate the open charm hadron abundances in e + e~~ 
collisions with an experimental input of the total open 
charm production. 

In this letter we consider the c.e. formulation |)| of 
the statistical coalescence model S. The experimental 
value of the J/tp multiplicity (J/4>) will be used to predict 
the open charm yield within the statistical coalescence 
model. We find a strong enhancement of the open charm 
production, by a factor of about 3, over the standard 
hard-collision model extrapolation from nucleon-nucleon 
to A+A collisions. 

In the grand canonical ensemble (g.c.e.) the thermal 
multiplicities of both open charm and charmonium states 
are given as (Bose and Fermi effects are negligible): 
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where V and T correspond to the volumqj and tempera- 
ture of HG system, rrij , dj denote particle masses and de- 
generacy factors and if 2 is the modified Bessel function. 
The particle chemical potential fj,j in Eq.(|l|) is defined as 



CjfJ-C 



(2) 



where bj , Sj , Cj denote the baryonic number strangeness 
and charm of particle j. The baryonic chemical poten- 
tial fiB regulates the baryonic density of the HG system 
whereas strange us and charm pc chemical potentials 
should be found from the requirement of zero value for 
the total strangeness and charm in the system (in our 



*To avoid further complications we use ideal HG formulae 
and neglect excluded volume corrections. 
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consideration we neglect small effects of a non-zero elec- 
trical chemical potential). 

In the c.e. formulation (i.e. when the requirement of 
zero " charm charge" of the HG is used in the exact form) 
the thermal charmonium multiplicities are still given by 
Eq.(|l]) as charmonium states have zero charm charge. 
The multiplicities (|l|) of open charm hadrons will, how- 
ever, be multiplied by the additional 'canonical suppres- 
sion' factor (see e.g. pp. This suppression factor is the 
same for all individual open single charm states. It leads 
to the total open charm multiplicity Nq in the c.e.: 



N c e = N 



h(Np) 
Io(N ) 



(3) 



where No is the total g.c.e. multiplicity of all open 
charm and anticharm mesons and (anti)baryons calcu- 
lated with Eq. (Q) and la , I\ are the modified Bessel func- 
tions. For large open charm multiplicity No >> 1 one 
finds h{No)/h{No) -» 1 and therefore N Q e -> N , i.e. 
the g.c.e. and c.e. results coincide. For No << 1 one has 
h(N o )/I (N o ) = N /2 and N Q e S- N ■ N /2, there- 
fore, Nq is strongly suppressed (like No/2, No << 1) 
in comparison to the g.c.e. result, No ■ 

Assuming the presence of the charm enhancement fac- 
tor 7 C the statistical coalescence model within the c.e. is 
formulated as: 
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when Nh is the total (thermal) multiplicity of parti- 
cles with hidden charm. Therefore, the baryonic num- 
ber, strangeness and electric charge of the HG system 
are treated in our approach according to the g.c.e. but 
charm is considered in the c.e. formulation where the 
exact charge conservation is imposed. 

We will proceed with Eq.(|4|) in the following way. As 
the (J/ip) multiplicities can be extracted from the NA50 
data on Pb+Pb collisions at 158 A-GeV for different val- 
ues of N p , we start from the requirement: 



7 C 



fjtot 
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to fix the 7 C factor. In Eq. 
multiplicity is calculated as 



the total J/ip thermal 



njtot 



Nj/^ + R(tf)N,y 
+R(Xi)N xl +R(X2)N 
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(6) 
(7) 



where N 



N^i, N Xl , N X2 are given by Eq.(|l|) and 



R(<ip') = 0.54, R(xi) = 0.27, Rfa) = 0.14 are the decay 
branching ratios of the excited charmonium states into 
J/ip. Eq.(H) will be used then to calculate the value of 
N^ r . This value will be considered as a prediction of the 
statistical coalescence model: the open charm yield has 
not yet been measured in Pb+Pb collisions at SPS. 



A reliable extraction of the J/tp yields from the pub- 
lished data appears to be non-trivialjj} The results for 
(J/ip) presented in Ref. || were evaluated from the data 
of the NA50 Collaboration |D| using the procedure de- 
scribed in p^j |. These results are presented in Table 1 to 
be used in the present analysis. 

We use the set of the chemical freeze-out parameters 



T = 168 MeV, hb = 266 MeV 



(8) 



They are fixed by fitting the HG model to the hadron 
yield data in Pb+Pb collisions at 158 A-GeV (the in- 
clusion of open charm and charmonium states does not 
modify the rest of the hadron yields). For a fixed num- 
ber of participants, N p , the volume V is calculated from 
N. p = Vub, where ng = Ub(T, hb) is the baryonic den- 
sity calculated in the g.c.e.. With the chemical freeze-out 
parameters (|J) we find the Nj"^ and No values using 
Eq.(0). Then we calculate 7 C from Eq.(|). Finally, Nfi r 
is calculated from Eq.(g). Note that T = 170 MeV leads 
to the HG value of the thermal ratio of (i/j') / (J/ip) — 0.04 
in agreement with data [|l6| at N p > 100. This fact was 
first noticed in Ref. JL7[ . 

The model results for central Pb+Pb interactions at 
158 A-GeV (N p = 100 ~ 400) are presented in Table 1. 
Assuming that N^ r scales as N£ we find a = 1.7. This 
value is larger than a = 4/3 expected in the hard-collision 
model. 

pQCD inspired models suggest values of N^ r — 0.15-^ 
0.3 in central Pb+Pb collisions at 158 A-GeV (a value of 
N^ r 0.17 is estimated in Ref. [g). The results pre- 
sented in Table 1 correspond to an enhancement of the 
open charm production by a factor of about 2-4. Al- 
though the values of Nj°fp, No and 7 C are rather sensitive 
to the temperature parameter, the model predictions for 
iVcg r remain essentially unchanged when T = 170 ± 10 
MeV is used. Note that a recent analysis of the dimuon 
spectrum measured in central Pb+Pb collisions at 158 
A-GeV by NA50 Collaboration Q suggests a significant 
enhancement of dilepton production in the intermediate 
mass region (1.5^-2.5 GeV) over the standard sources. 
The primaryn interpretation attributes this observation 
to the enhanced production of open charm Jig] : about 
3 times above the pQCD prediction for the open charm 
yield in Pb+Pb collisions at SPS. This value is in agree- 
ment with our model analysis. Moreover, the analysis of 
the centrality dependence of the dimuon excess was done 
in Ref. 0] assuming that additional dimuon pairs come 



'We are thankful to M. Gazdzicki and K. Redlich for useful 
comments. 



Alternative explanations are also suggested jl9| 53. 
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entirely from decays of D and D. The open charm yield 
was found to be proportional to N p 7 . This precisely co- 
incides with our result. 

A possible source of the open charm enhancement in 
A+A collisions with respect to the direct extrapolation 
from p+p data may be the broadening of the phase space 
available for the open charm due to the presence of the 
hadron and/or the quark-gluon medium: 

Let a cc pair with an invariant mass below the open 
charm threshold, 2Mb = 3.7 GeV, be created in a p+p 
collision. In the vacuum it cannot hadronize to an open 
charm hadron pair and thus it has to be transformed into 
non-charmed states to respect energy-momentum conser- 
vation. In A+A collisions, however, such a pair has a 
chance to shift its invariant mass above the open charm 
threshold due to secondary interactions with hadron 
and/or quark-gluon medium, thus resulting in additional 
open charm production. The strongest enhancement of 
the open charm production should be expected in the 
quark-gluon plasma. Due to the Debye screening, c and 
c behave like free particles. Therefore, almost all sub- 
threshold cc-pairs created in the hard collisions hadronize 
into open charm mesons and baryons at the later stage 
of the A+A reaction. This effect can be substantial at 
moderate center-of-mass energies per nucleon, while at 
higher energies the threshold effects becomes less impor- 
tant. A quantitative estimate of this effect is under way 
PH . The above mechanism does not alter the production 
of direct Drell-Yan pairs: both the threshold effects and 
the interaction with the nuclear medium are negligible in 
that case. 

In conclusion, the statistical coalescence model with 
exact charm conservation has been formulated. The 
canonical ensemble suppression effects are important for 
the thermal open charm yield even at N p = 400. These 
effects become crucial when the number of participants 
N p decreases. From the J/ip multiplicity data in Pb+Pb 
collisions at 158 A-GeV the open charm yield is predicted: 



N ai 



dir 



0.6 -7- 0.7 



in central collisions. This is about a factor 3 above the 
pQCD prediction for the open charm yield in Pb+Pb 
collisions at SPS. The statistical coalescence model pre- 
dicts also an N p dependence of N^ r ~ Np 7 , which is 

stronger than the standard dependence, Np^ 3 , in the 
hard-collision model. These predictions of the statisti- 
cal coalescence model can be tested in the near future at 
the CERN SPS, where measurements of the open charm 
are planned. 
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2.2 ± 0.2 
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